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ITERATIVE  CURVE  FIT  OF  A  POWER  LAW  EGR.  ATION  TO 
SOLID  PROPELLANT  SHIFT  FACTOR  D.vTA 


1.  INTRODUCTION 

The  JANNAF  Solid  Propellant  Structural  Integrity  Handbook  (Ref.  I)  suggests  an 
empirical  power  law  relationship  between  temperature  (T)  and  the  time-temperature 
shift  factor  (ay).  To  supplement  existing  solid  propellant  thermal  transient  viscoelastic 
response  software  (Ref.  2)»  an  iterative  method  of  fitting  th ;  transcendental  equation 
relating  T  and  ay  to  a  given  set  of  experimental  data  was  de /eloped.  This  method  was 
implemented  in  software  formulated  for  use  on  a  Hewlett-Pac<ard  981 5A  programmable 
calculator  (HP  98I5A).  It  was  decided  that  with  minor  charges  the  software  could  be 
modified  to  fit  an  arbitrary,  single-valued,  algebraic  or  transcendental  equation  to 
experimental  data  with  the  number  of  unknown  constants  li  nited  only  by  the  storage 
capacity  of  the  calculator. 

2.  SUMMARY 

The  computer  program  developed  for  the  HP  98I5A  successfully  determines  values 
for  the  unknown  constants  of  the  power  law  ay  equation  in  fen  cycles  or  less,  with  a 
total  run  time  of  approximately  ten  minutes.  These  constants  >roduce  a  curve  which  fits 
a  given  set  of  log  ay  vs  T  data  with  less  error  than  the  trial  and  error  method  suggested 
by  the  JANNAF  Solid  Propellant  Structural  Integrity  Handbook  (Ref.  I).  It  also  produced 
a  better  fit  than  the  direct  approximate  solution  suggested  in  A  FRPL-TR-8 1  -80  (Ref.  2). 

This  program  was  also  expanded  to  fit  a  general  function  (either  algebraic  or 
transcendental,  and  having  no  more  than  three  undetermined  c  instants)  to  a  given  set  of 
data.  Due  to  lack  of  knowledge  in  making  initial  guesses  for  tl  e  constants,  this  curve  fit 
does  not  converge  as  rapidly  as  the  specific  curve  fit  for  the  power  law  ay  equation.  In 
fact,  if  the  constants  chosen  produce  too  large  an  error,  the  iterative  process  will  fail 
and  new  trial  values  must  be  selected.  Modification  of  the  storage  capacity  of  this 
program  would  allow  equations  with  more  than  three  unknown  c  instants  to  be  evaluated. 


3.  SIGNIFICANCE  OF  A  POWER  LAW  REPRESENTATION  -  W  LOG  (aT)  VS  T  PLOT 


3.1  Definition  of  (crrj 

When  a  viscoelastic  tensile  test  specimen  is  strelshed  (at  a  constant  speed 
under  uniform  ambient  temperature)  to  a  specific  strain  level,  the  force  required  to 
maintain  that  strain  level  thereafter  decreases  with  time.  Ttv*  modulus  of  elasticity  will 
also  be  time  dependent,  and  is  denoted  as  the  relaxation  rrodulus  E(t).  This  test  is 
referred  to  as  a  relaxation  test,  and  if  it  is  performed  at  a  seri  ?s  of  ambient  temperature 
levels,  a  group  of  similar  shaped  modulus  curves  plotted  as  log  E(t)  vsiog  time  (t)  may  be 
obtained. 


The  assumption  is  now  made  that  solid  propellants  are  thermorheo logically 
simple  materials  (Ref.  3)  and  that  the  time-temperature  superposition  principle  is 
applicable.  This  implies  that  the  group  of  curves  generated  fr<  m  the  relaxation  tests  are 
portions  of  a  single  master  relaxation  curve  displaced  (shifted  along  the  log  (t)  axis  due 
to  the  different  temperature  levels  at  which  the  tests  were  run. 

The  time-temperature  shift  factor  (aj)  is  defir  ed  as  a  non-dimensional 
measure  of  the  shift  in  the  log  E(t)  vs  log  (t)  curve  due  to  a  c  xmge  in  the  temperature 
of  the  relaxation  test.  Quantitatively,  aT  is  the  distance  along  the  log  (t)  axis  between  a 
value  of  log  E(t)  an  a  reference  modulus  curve  and  an  equal  value  of  log  E(t)  on  the 
modulus  curve  for  the  temperature  of  interest.  In  practice,  the  relaxation  modulus 
curves  are  treated  as  a  single  "master  curve"  of  E(t)  vs  log  %  where  £  t/ay  is  defined 
as  the  "reduced  time"  variable.' 

3.2  Power  Low  qt  Equation 

An  empirical  power  law  equation  has  been  develop  ?d  (Ref.  2)  to  determine  a 
value  of  aj  for  a  relaxation  test  run  at  a  temperature  abo  'e  or  below  an  arbitrary 
reference  temperature  (Tr)  with  the  restriction  that  the  ten  perature  of  interest  must 
be  greater  than  the  glass  transition  temperature  of  the  propel  It  nt.  The  relationship  is  as 
follows: 
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Ta  is  a  positive  or  negative  constant  with  a  suggested  value  of  I0°F  to  20°F 
below  the  gloss  transition  temperature,  and  m  is  a  positive  constant  (Ref.  2).  The  value 
of  Tr  is  usually  in  the  range  of  70°F  to  80°F  (Ref.  2). 

3.3  Application  of  the  Power  Law  ot  Equation  to  Struc  oral  Analysis 

The  reduction  of  a  thermal  transient  problem  to  an  "ordinary  linear 
viscoelastic"  problem  (Ref.  3)  is  based  on  the  Moreland  Lee  hypothesis  used  in 
conjunction  with  the  previous  assumption  that  the  materi  il  is  thermorheologically 
simple.  The  Moreland-Lee  hypothesis  asserts  that  the  reduced  time  under  variable 
temperature  conditions  is  given  by 


The  power  law  representation  of  ay(T)  was  developed  to  provide  an 
analytically  integrable  function  for  the  case  of  a  linear  temperature-time  history  and  a 
restricted  class  of  nonlinear  temperature-time  histories. 

4.  PROBLEM  STATEMENT 

From  the  preceding  section,  it  becomes  evident  that  the  analyst  needs  a  method  to 
determine  the  values  of  the  arbitrary  constants  Ta  and  m  in  the  power  law  ay  equation  in 
order  to  apply  the  Moreland-Lee  hypothesis.  Data  relating  <r  and  T  may  be  generated 
from  a  group  of  relaxation  tests  performed  over  a  range  of  temperatures.  The  least 
squares  method  of  fitting  a  curve  to  determine  the  values  <>f  the  unknown  constants 
results  in  a  set  of  normal  equations  non-linear  in  the  constant;  Ta  and  m.  Therefore,  a 
method  is  needed  to  generate  the  values  of  Ta  and  m  for  af  vs  T  data  such  that  the 
power  law  ay  equation  fits  this  data  with  minimum  error;  i.< , to  fit  a  transcendental 
equation  with  two  or  more  unknown  constants  to  a  given  set  of  lata. 


The  term  "error"  will  have  two  definitions  in  the  fotio'/ing  discussion.  The  root 

A 

mean  square  error,  E,  is  defined  as: 

l 


1-1 


(3) 


where  n  is  the  number  of  data  points  and  AYj  Is  the  distance  in  the  y  direction  between  a 
data  point  (Xi,  Yj)  and  the  point  on  the  curve  corresponding  to  Cj. 

The  sum  of  the  squared  errors,  E,  is  defined  as: 

E  *  (‘  V2  «) 


5.  SOLUTION 

5.1  Assumptions 

The  following  method  of  transcendental  curve  f  tting  was  based  on  two 
assumptions: 

(«)  The  n  unknown  constants  C;,  i  =  l...n  produce  an  (n  +  I)  dimensional  error 
surface  which  represents  the  error,  E,  in  the  curve  fit  du.;  to  a  particular  set  of 
constants.  When  the  intersection  of  level  planes  with  this  sur  ace  are  projected  into  rv 
space,  they  are  assumed  to  produce  a  group  of  inscribed  contour  lines  (see  Figure  I). 
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Figure  1.  Hypothetical  Error  Surface,  E  *  ftCj),  for  Two  Unknowi  Constants. 


(2)  The  function  describing  the  error  surface,  E  =  f(C|),  with  respect  to  any 
constant  C;  may  be  accurately  approximated  over  a  small  interval  by  a  quadratic 
equation,  f(C;)  =  A(Cj)2  +  B(Cj)  +  D,  and  a  minimum  value  for  his  quadratic  exists  in  the 
neighborhood  of  that  small  interval.  The  plots  in  Figures  7  and  3  are  presented  as 
examples  of  Hie  curve  representing  E  =  f(C{)  developed  by  fitting  a  quadratic  through  the 
three  points  marked  "+."  The  points  marked  "o"  are  actual  er  ror  points;  it  can  be  seen 
that  over  the  interval  used  to  develop  the  curve,  the  quadratic  accurately  approximates 
the  actual  error  curve  and  that  a  minimum  value  exists  in  or  ne  ir  this  interval. 
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The  data  for  Figures  2  and  3  were  generated  fro.;*  the  log  ay  vs  T  data  used 
in  the  sample  problem  of  Appendix  A. 

If  the  trial  values  of  the  constants  are  not  relative  /  close  to  the  values  that 
minimize  the  error  in  the  fit,  the  second  assumption  will  not  be  valid.  Also,  Figure  3 
shows  that  points  on  the  actual  error  curve  not  in  the  imme  Hate  neighborhood  of  the 
points  used  to  generate  the  quadratic  are  not  necessarily  well  p  edicted  by  the  quadratic. 

5.2  Selection  of  initial  Trial  Values 

The  curve  fit  technique  developed  is  an  iterative  pr  ocess  that  requires  initial 
guesses  at  the  unknown  constants.  The  technique  is  very  «nsitive  to  these  initial 
guesses.  Therefore,  a  method  was  devised  to  develop  initial  gjesses  for  Ta  and  m  from 
ay  vs  T  data. 

If  a  value  of  Ta  equal  to  Tg  is  substituted  into  equr  tion  (I),  a  value  of  m  may 
be  determined  from  a  least  squares  fit  of  the  linear  relationship  (see  Figure  4): 


tt  T*  -  T«1  /  t  T  - 


Figure  «.  Least  Squares  Fit  of  Log  vs  LoqpR  ~  V[  for  staple  i  "Oblea  In  Appendix  A. 


A  user  knowing  the  values  of  Tr  and  Tg  can  thus  obtain  fairly  occurate 
guesses  at  Tq  and  m.  If  Tr  is  unknown,  a  value  between  70of  and  80°F  should  be  used 
based  on  the  definition  of  Tr  in  Section  3.2.  If  Tg  is  ur  known,  a  trial  and  error 
procedure  of  guessing  values  for  Tg  and  plotting 


should  be  used  until  a  nearly  straight  line  plot  is  obtained.  Tne  guess  of  Tq  that  yields 
this  straight  line  plot  should  then  be  used  as  the  initial  guess  fo*  Ta. 

The  JANNAF  Solid  Propellant  Structural  Integrity  Handbook  (Ref.  I)  suggests 
the  above-described  trial  and  error  procedure  to  obtain  the  fin  il  ax  equation  assuming  a 
value  of  Ta  1 0°F  to  20°F  below  Tg  (Ref.  3).  However,  it  was  found  in  the  present  study 
that  the  best  fits  were  obtained  using  the  above  method  to  <  enerate  guesses  only  and 
then  iteratively  operating  on  these  guesses  to  obtain  values  t  >at  minimize  the  error  of 
the  fit. 


5.3  Iterative  Optimization  of  the  Constants 

The  initial  guesses  of  the  unknown  constants  descri  >e  a  point  A  (see  Figure  5) 
in  the  Ta  -  m  plane.  The  value  of  m  is  temporarily  held  constant,  and  points  B  and  C  are 
generated  by  adding  and  subtracting  a  user -specified  multiple  p)  of  the  present  value  of 
Ta  to  or  from  Ta,  respectively.  This  search  interval  is  labeled  ,\Ta,  and  is  equal  to  C.  = 
p  C . ,  where  C;  equals  Ta  in  Figure  5.  In  subsequent  cycles,  t'lis  search  interval  will  be 
multiplied  by  a  decreasing  factor  generated  from  the  function 


(1  -  6)  +  ee(1'b)a 


(7) 


V  '  »  ITT 


where  b  is  the  number  of  the  current  iteration  cycle,  and  o  and  $  are  user  specified 
constants  that  determine  the  type  (linear,  exponential,  et-:.)  and  rate  of  interval 
reduction.  This  reduction  in  the  search  interval  allows  more  accurate  fits  of  a  quadratic 
approximation  to  the  error  curve,  producing  a  more  accurate  approximation  of  the 
minimizing  value  of  C;  as  the  iterative  process  gets  closer  to  the  actual  minimum  values. 
In  a  more  general  sense,  then,  AC;  may  be  defined  as 

ac,.  =  (P)(c.)  jji-e)  +  ee(1'b)a]  (8) 


To  assure  that  the  quadratic  approximation  to  the  error  has  a  minimum,  the 
search  interval  will  be  expanded  first  in  the  negative  C;  direction  and  then  in  the 
positive  C;  direction  until  the  values  of  E  associated  with  the  >oints  B  and  C  (E(B),  E(C)) 
are  greater  than  E(A).  In  the  example  shown,  Ta  is  expended  in  the  negative  Tn 
direction  until  the  point  C  is  located  such  that  E(C)  >  E(A). 


A  quadratic  is  then  fitted  through  C,  A,  and  3,  and  the  first  partial 
3E(r  ,  m) 

derivative,  - 1 - ,  is  set  equal  to  zero  to  find  a  new  v  ilue  of  Ta  (point  D)  which 

<7  I  _ 


minimizes  the  error.  See  Appendix  B  for  the  methods  used  to  :it  the  quadratic  and  take 
the  partial  derivative  of  the  error. 


Since  the  quadratic  fit  is  only  an  approximation  to  fhe  actual  error  function, 
a  value  of  E  corresponding  to  point  D,  E(D),  is  calculated  and  compared  with  the  starting 
points  E(A).  If  E(D)  is  greater  than  E(A),  the  new  value  T'a  will  be  rejected.  If  E(D)  is 
less  than  or  equal  to  E(A),  the  value  T'a  will  replace  Ta  as  shown  in  Figure  5. 
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F1*ur»  $.  I  ten  tin  Optimization  of  tee  Constants,  Cycla  1. 


An  identical  process  now  operates  on  m  (Ta  ten  porarily  fixed),  with  the 
quadratic  fitted  through  points  D,  F,  and  G.  A  minimum  valui  of  the  error  is  found  at 
point  H  (T'q,  m').  The  value  of  E  corresponding  to  point  H  is  r  ompared  to  E(D),  and  the 
same  acceptance  criterion  for  this  new  point  is  applied.  In  Figure  6,  E(H)  was  less  than 
E(D). 

One  cycle  of  the  procedure  is  now  complete;  and  due  to  the  comparisons,  first  of 
E(A)  and  E(D)  and  then  E(H)  and  E(D),  it  is  evident  that  the  i  rocess  has  moved  in  the 
proper  direction  to  reduce  the  error  in  the  curve  fit. 

A  likely  starting  point  for  the  second  cycle  mfcht  be  point  H  since  it 
corresponds  to  the  minimum  E  at  present.  However,  to  avcid  the  possibility  of  the 
process  becoming  "stuck"  in  a  local  valley  of  the  error  surf  ice,  a  point  A',  midway 
between  the  newly  calculated  minimum  point  H  and  the  staiting  point  of  that  cycle 


(j-oint  A),  is  determined  and  the  process  is  started  over  (see  Figure  6)  at  this 
intermediate  point.  The  next  cycle  produces  a  new  minimum  at  point  H'  and  a  new 
intermediate  point  A".  Since  the  minimum  value  of  E  and  the  corresponding  values  of  Ta 
aid  m  are  always  stored  in  memory,  these  values  can  be  (and  are)  selected  on  the  final 
c  /cle  rather  than  using  the  constants  associated  with  A'  for  the  final  values. 


n*«r«  t.  Itaratlvi  0»t  Mutton  of  th*  ComUnts,  Cycle  2. 


This  process  is  repeated  for  a  user-specified  numtor  of  cycles  or  until  the 
value  of  E  at  any  point  in  the  process  becomes  less  than  a  user-  specified  minimum  E. 

It  should  be  noted  that  Tr  is  considered  a  known  constant,  and  the  power-law 
a**  equation  should  ideally  pass  through  the  point  (Tr,  0)  on  a  If  g  ay  vs  T  plot  since 


However,  the  log  ay  vs  T  data  will,  most  likely,  lw  e  inherent  error,  and  the 
ctrve  which  best  fits  this  data  may  not  necessarily  pass  through  (Tr,  0),  requiring  the 
procedure  to  iterate  on  a  third  constant,  Tr.  Since  both  the  assumptions  of  Section  5.1 
held  for  an  n-dimensional  problem,  and  since  the  equation  fir  distance  between  two 
points  may  be  extended  to  n  dimensional  space,  this  third  cor*  tant  presents  no  problem 
even  though  it  produces  a  four  dimensional  error  surface. 
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Appendix  A  illustrates  a  sample  problem  that  uses  this  iterative  process. 
Figure  A- 1  shows  a  plot  of  the  curve  generated  by  the  iterative  process  along  with  the 
data  used  in  the  problem. 

6.  GENERAL  CURVE  FITTING 

The  method  of  optimizing  the  constants  of  a  function  de?  bribed  in  Section  5  may  be 
applied  to  any  single-valued  equation,  either  algebraic  or  trot  scendental.  The  problem 
again  arises  of  picking  initial  values  for  the  unknown  conston  s.  If  arbitrary  values  are 
selected,  the  user  must  take  care  in  choosing  those  values  •  uch  that  the  equation  of 
interest  can  be  evaluated  over  a  specified  range  when  calculating  the  error.  Also,  the 
program  must  be  modified  to  allow  for  the  possfcle  nonexistence  of  a  minimum  during 
any  particular  variable  search. 

Three  modifications  were  made  to  the  program  previr  jsly  developed  to  fit  the 
power  law  aj  equation.  First,  the  portion  of  the  program  that  loads  the  data  was 
modified  to  receive  guesses  rather  than  generate  initial  valu  »  from  the  data  entered. 
Second,  the  user  must  program  the  equation  of  interest  nto  the  subroutine  that 
calculates  the  sum  of  the  squared  errors.  Third,  the  portion  jf  the  program  that  fits  a 
quadratic  through  the  error  values  and  minimizes  the  function  was  modified  to  take  into 
account  the  possibility  that  the  curve  has  no  minimum.  I  the  curve  is  seeking  a 
maximum,  the  program  defaults  to  the  one  of  three  poinl  ?  used  to  determine  the 
quadratic  with  the  smallest  E  value  and  then  proceeds  to  iterat  •  on  the  next  constant.  If 
the  curve  is  seeking  a  minimum  outside  the  interval  (C;-nCj,  C  ;  ♦  nCj),  the  interval  over 
which  Cj  is  expanded  is  doubled  both  in  the  plus  and  min's  Cj  direction  until  the 
minimum  value  of  Cj  lies  within  the  interval  of  (Cj  -  nCj,  Cj  ►  nCj)  with  the  increased 
value  of  n. 

Since  arbitrary  initial  guesses  are  being  used,  the  conven  ence  of  the  general  curve 
fit  is  not  as  rapid  as  the  power  law  aj  fit.  However,  successi\  3  use  of  the  program  will 
allow  a  user  to  narrow  down  the  choices  for  the  constants  or  d  finally  arrive  at  values 
that  are  accurate  enough  for  assumption  (2)  of  Section  5.1  to  ac  jin  be  valid. 
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APPENDIX  A.  SAMPLE  PROBLEM 


This  Appendix  demonstrates  a  sample  curve  fit  on  the  HP-981 5A  for  the  aj  vs  T 
data  tabulated  at  the  end  of  the  Appendix.  Statements  in  quotations  are  user  cues 
printed  by  the  program.  Numbers  and  words  shown  in  boxes  correspond  to  keystrokes 
needed  to  run  the  program. 

I.  Subroutine  to  Load  Data 

A.  User  enters  the  following  to  start  the  program: 

1  jEnterl  fT|  {Load  I  IcT  ear  ]  fEnd]  fR7sl 

R.  "PROGRAM  TO  LOAD  DATA  FOR  A(T)  VS  T  CUR\  E  FIT” 

"NO.  OF  DATA  PTS?" 

UiiZsl 

(User  enters  the  number  of  data  points  (up  to  twent  ').  If  a  number  greater 
than  twenty  is  entered,  the  message  will  be  repeate  I.) 

C.  ’ENTER  DATA  (X,Y)” 

Enter  coordinates  of  data  points  (X  is  T  in  °F;  Y  is  log  aj). 
f-45{  |R/Sj(actual  by  stroke  sequence  is  [3J  fCHSl  jR/SQ 
"-45" 

e&i  m 

"7.52" 


n52llR7s} 

"152” 

E MM] 

"-2.30" 

D.  "ENTER  REF.  TEMP.  DEGREES  F" 

(If  this  value  is  unknown,  see  Section  5.2  for  a  range  of  values.) 


E.  "ENTER  GLASS  TRANSITION  TEMP.,  DEGREES  F" 


(This  is  used  as  an  initial  guess  for  Ta.  If  Tg  is  unknown,  see  Section  5.2  for  a 
method  of  generating  this  value.) 

mm 

"-92" 

A  trial  value  of  m  is  now  computed  by  a  least  squares  fit  of  log  ay  vs 
log  (Tr  -  Tq)  data,  m  being  the  slope  of  this  line: 

"m  =  13.087" 

F.  "STORE  ON  TAPE?" 

|R/51(For  yes;  |tiny  numeric  keyj  then  jR/S}  for  no.) 

This  records  the  data  points  as  well  as  the  constants  on  tape. 

G.  "LOAD  CURVE  FIT?" 

IR/Sl  (For  yes:lany  numeric  key^  then  IR/Sl  for  no.) 

This  loads  the  subroutine  that  iteratively  optimizes  the  constants. 


Iterative  Optimization  of  Constants 


A.  If  this  subroutine  has  not  been  loaded  by  the  subroutine  to  load  data,  enter: 
j0|  | enter)  f?)  jloadl  ICIeoH  lEndl  fR/Sl 
(otherwise  begin  with  Step  B.) 


3. 


"NO  OF  CONSTANTS  3  MAX" 


Enter  number  of  constants  in  the  equation. 

mm 


c. 


•ENTER  SEARCH  INTERVAL" 

mm  /Tr  -  T\ 

(Since  ay  =  m  loq  — j-  —  y  )  must  be  evaluated  at  points  generated  by  this 
constant  (p),  care  must  be  taken  in  selecting  this  constant  so  the  equation 
may  be  evaluated  over  the  interval.  In  this  problem,  p  is  limited  to  .489.  A 
value  of  p  =  .489  would  cause  the  program  to  take  the  log  of  a  negative 
number  when  iterating  on  Tai  producing  an  error  condition.  A  value  of  0.1  or 
smaller  is  recommended.) 


D.  TiNTER  B" 

mm 

E.  "ENTER  ALPHA" 

MM 

(These  values,  ALPHA  and  BETA  (B),  control  the  reduction  of  the  interval 
with  each  cycle.  B  =  .9,  ALPHA  =  .2  give  an  exponential  decay  to  a  constant 
value.) 

F.  "ENTER  TOLERANCE" 

i mm 

(This  value  corresponds  to  a  minimum  E  that  the  iterative  process  is  trying  to 
achieve.  If  this  value  is  reached,  the  program  will  stop  and  print  the 
constants  associated  with  it.) 

G.  "NO  OF  ITERATIONS" 

Pol  IrTsI 

(This  corresponds  to  the  number  of  cycles  to  be  nm.  Convergence  to  three 
decimal  places  is  usuolly  achieved  in  10  cycles  or  l-;ss  with  each  cycle  taking 
approximately  15  seconds  if  all  three  constants  are  iterated  on.) 

H.  "SET  FLAG  I  TO  CONSTRAIN  C|" 
fR/S? 

Setting  Constrain  s 
Flag  I  Ta 

Flag  2  Tr 

Flag  3  m 

The  best  fit  is  obtained  by  leaving  all  three  constants  unconstrained.  If  Tr 
was  to  be  constrained,  the  user  would  press  ISFGl  Q  j  jR/S} 

I.  "INITIAL  VALLE  OF  CONSTANTS" 

"92"  (NOTE:  The  program  assumes  Ta  is  negative  and  only  uses  the 
magnitude  of  this  number.) 


BE 


"13.067" 


Each  constant  is  now  redisplayed  on  the  screen.  If  that  constant  is  to  be 
used,  press  [R/^.  If  a  different  value  is  to  be  used,  enter  the  value  Cj  and 
press  lR/5t 

NOTE:  When  the  final  constant  is  printed,  and  (R/S)  is  hit,  the  |SFG|  key  should  be  hit 
immediately  if  the  user  wishes  to  have  the  error  values  printed  out.  Re-hitting jSFCO will 
s  op  the  values  from  being  printed. 


J.  Output  and  User  Termination 

The  value  of  E  for  each  cycle  is  printed,  followed  by  the  value  of  E  after  an 
iteration  on  one  of  the  constants  is  performed,  assuming  fSFG)  is  pressed  in  Step  I.  The 
program  will  run  the  specified  number  of  cycles  unless  E  at  any  point  becomes  less  than 
the  tolerance,  or  the  user  terminates  the  program  by  pressing |R/S| Q.  If  the  program 
runs  its  full  number  of  cycles,  the  values  of  the  final  E  and  corresponding  constants,  as 
well  as  the  minimum  E  and  the  corresponding  constants,  will  be  printed.  If  the  user 
terminates  the  program  at  some  intermediate  point,  the  number  of  the  current  iteration 
cycle  is  printed  along  with  current  E  and  constants  as  well  as  minimum  E  and  constants. 
Finally,  if  the  tolerance  is  met,  the  value  of  E  and  constants  which  met  or  exceeded  the 
tolerance  will  be  printed.  The  number  which  flashes  on  the  display  indicates  the  current 
iteration  cycle.  Figure  A-l  shows  a  plot  of  the  data  used  for  the  problem  and  the  curve 
fit  obtained  by  the  iterative  method.  The  input  data  is  also  tabulated  below. 


Data  used  in  sample  problem 


Log  aj 

T  (°F) 

7.52 

-45 

6.60 

-42 

5.19 

-29 

3.84 

-2 

2.20 

19 

0.00 

76 

-1.23 

i  18 

-2.30 

152 
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APPENDIX  B.  CURVE  FIT  AND  MINIMIZATION  OF  ERROR 


The  following  discussion  outlines  the  procedure  used  to  simultaneously  fit  a 
quadratic  through  the  three  points  used  to  generate  the  error  •  urve  and  find  a  value  of  C 
that  minimizes  this  function  (see  Figure  B-l).  In  this  discussi  >n,  only  one  of  the  Ci*s  is 
considered  to  vary. 

H  =  A(C|)2  +  B(C|)  +  D  =  E(C|)  (B-l) 

I  =  A(C2)2  +  B(C2)  +  D  =  E(C2>  (B-2) 

J  =  A(C3)2  +  B(C3)  +  D  =  E(C3)  (B-3) 

Simultaneous  solution  of  equations  (B-i)  through  (B-3)  yi  -Ids  a  quadratic  fit  of  the 
squared  error  E(C)  =  E  ( A  Y;)2  through  the  three  points. 

Cj  (J  -  1)  +  C2  (H  -  J)  +  C3  (I  -  H)  (B-4) 

A  -  9.9  O  6  "  1  _____ 

C1  (C3  -  C2 >  +  C2  (C1  *  C3  >  +  C3  ^C2  -  ci  ) 


B  = 


(I  -  J)  -  A(C„2  -  C,2) 

- 


(B-5) 


The  function  A(Cj)2  +  B(C;)  +  D  describes  the  squared  ern  r,  E(Cj).  Taking  the  first 
partial  derivative  of  E  with  respect  to  C;  and  setting  it  eqix  I  to  zero  determines  the 
minimum  value  of  C| 


+  BCi 

IT 

l 


♦  0) 


0 
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0  =  2AC{  +  B 

an^  ^-i  min  =  — 
2A 


(B-6) 


Direct  substitution  of  equations  (B-4)  and  (B-5)  into  equation  (B-6)  yields  the  value 
of  C;  that  minimizes  the  squared  error  E.  Due  to  the  location  of  the  three  points  C|, 
C2,  and  C3,  the  value  computed  in  equation  (B-6)  is  assured  to  be  a  minimum.  See 
Section  5.3  for  an  explanation  of  the  process  of  locating  C|,  C2,  and  C3. 


£  *  f(C,)  - 
A{C,)Z+B(Ci)*D 


Flgwrt  8-1.  Ctirvt  Fit  tmi  Rlufalzotlen  of  £  ■  f(C.) 
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